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Thirty  three  crack  growth  curves  were  analysed  in  order  to  estimate 
the  logarithmic  variance  of  fatigue  crack  growth  rate  for  several  7XXX 
Aluminium  Alloys  in  various  conditions.  The  estimated  Standard  Deviation 
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Thirty  three  crack  growth  curves  were  analysed  in  order  to 
estimate  the  logarithmic  variance  of  fatigue  crack  growth  rate  for  several 
7 XXX  Aluminium  Alloys  in  various  conditions.  The  estimated  Standard 
Deviation  common  to  all  of  these  is;  s  *  0.1241  with  18  degrees  of  freedom 
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1. 


INTRODUCTION 


In  most  reliability  analyses  for  fatigue  life  prediction 
the  form  of  crack  growth  is  deterministic  in  that  there  is  no  random¬ 
ness  in  the  assumed  length  of  any  crack  that  has  grown  for  a  given 
number  of  cycles.  This  allowance!  now  has  been  made  in  recent  work  and 
the  associated  computer  programmes* 

However  there  is  little  convenient  data  about  rate 
scatter  known  to  the  author  in  the  literature.  Although  there  is  data 
that  can  be  used  it  is  inevitably  in  graphical  form  which  would  entail 
great  effort  for  any  statistical  analysis.  Some  of  this,  referring  to 
d6ac  steel,  has  been  analysed  by  Jost  and  Esson2. 

For  the  7XXX  series  of  aluminium  alloys  the  author  has 
recently  studied  an  Alcoa  report3  containing  unequivocal  data  of  a 
manageable  form  and  size  (33  specimens)  covering  a  wide  range  of 
conditions.  Although  the  amount  of  useable  data  is  still  small,  present 
circumstances  justify  its  analysis  for  scatter. 

2.  DATA 


This  consists  of  plotted  crack  growth  between  1.0  and 
3.5  cm  in  flat  tension  specimens  102  mm  x  6.4  tm  in  cross  section.  For 
each  specimen  two  crack  grov'th  curves  were  supplied  for  growth  from 
either  side  of  a  slit.  It  was  found  that  these  usually  coincided. 

This  is  an  important  result  in  itself  because  it  indicates  that  scatter 
in  crack  growth  arises  mostly  from  differences  between  snecimens  and 
not  during  the  growth  of  a  particular  crack. 

Replication  among  this  data  is  never  more  than  4  and 
usually  2.  The  main  comparisons,  apart  from  microstructure  are  between 
7010,  7050,  7075  and  7475  alloys  in  T6  and  T7  tempers.  These  are 
confounded  with  differing  load  programmes,  the  majority  being  constant 
amplitude  with  added  overload  cycles.  Table  1  shows  the  data  finally 
used  together  with  Figure  numbers  from  Reference  3  and  distinctive 
features  of  each  test.  The  unit  of  life,  centimetres,  is  explained 
below.  The  results  from  Figures  A12  and  A30  are  taken  at  a  crack  length 
of  2.28  cm. 

3.  ANALYSIS 


The  growth  rate  of  any  crack  depends  upon  local  stress 
intensity  and  other  influences  which  may  be  grouped  under  the  heading 
of  material  parameters.  The  coincidence  of  crack  growth  curves  mentioned 
above  indicates  that  the  latter  are  constants  for  any  particular 
specimen. 


Perusal  of  most  data  for  da/dN  vs  AK  on  logarithmic  scales 
indicates  that  variance  of  log (da/dN)  is  sensibly  constant  while  its 
average  follows  the  well  known  sigmoidal  curve.  Since  independence  of 
mean  and  variance  is  best  characterised  by  the  Normal  distribution  this 
suggests  that  log (da/dN)  is  Normally  distributed;  to  assess  airworthi¬ 
ness  one  should  estimate  this  logarithmic  variance.  In  these  conditions 
it  can  be  proved  that  the  number  of  cycles  or  programme  blocks  applied 
during  the  time  a  crack  grows  between  two  lengths  is  also  log-Normally 
distributed  with  the  same  variance.  Accordingly  for  each  specimen  the 
cycles  or  number  of  programmes  occupied  in  crack  growth  from  1.0  cm  to 
3.5  cm  were  measured  from  the  relevant  figure  by  a  centimetre  rule. 

Some  typical  figures  and  a  measurement  are  reproduced  here.  Because  the 
graphs  were  standardised  and  variation  of  log-times  were  the  only 
quantities  of  interest  the  measured  distance  could  be  used  directly. 

This  was  allowed  by  the  fact  that  different  scales  corresponded  to 
significant  differences3  caused  by  alloy  specification ,  tender  or  load 
programne.  In  terms  of  analysis  of  variance  the  estimators  of  error 
variance  are  not  affected  even  though  calculated  block  or  treatment 
effects  may  become  quite  meaningless.  Differences  due  to  deliberate 
changes  in  microstructure  were  taken  as  errors  if  the  alloy- temoer- 
load  designation  remained  the  same. 

3.1  Computations 

Data  from  Table  I  was  first  grouped,  as  in  Table  II,  into 
the  experimental  designs  shown.  Analyses  of  variance*  on 
logarithms  of  each  of  these  then  supplied  four  error  sums  of 
squares.  Being  conpatible,  these  were  pooled.  The  significance 
of  treatment  and  block  effects  confirmed  the  results  of  Reference 
3.  In  Table  II  the  treatments  and  blocks  are  indicated  by 
referring  to  the  first  entries  from  Table  I. 

The  appendix  illustrates  the  analysis  of  the  non- 
orthogonal  design  from  the  first  two-way  classification  in  Table 
II.  Crout9s  method  has  been  used,  adapted  to  a  pseudo  "inversion 
to  cope  with  the  reduced  rank  of  the  model. 


4.  CONCLUSIONS 

The  estimated  logarithmic  standard  deviation  from  the 
useable  data  in  Reference  3  is;  s  *  0.1241  with  18  degrees  of  freedom. 
This  estimate  holds  for  differing  types  of  7 XXX  aluminium  alloys* 
constant  amplitude  and  block  programme  loading.  The  errors  it 
swmarises  include  effects  of  microstructure*  heat  treatment  and  purity 
not  described  by  the  alloy  specification. 
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TABLE  I 


DATA  USED  FROM  REFERENCE  3 


Figure 

Life  (cm) 

Alloy 

Load  Type 

Remarks 

A7 

9.9 

7010-T7 

CA 

hi  purity 

A8 

10.5 

7475-T7 

CA 

A9 

11.2 

7475-T7 

CA 

HI  CU 

A10 

8.5 

7010-T7 

CA 

All 

12.8 

7050-T7 

CA 

LO  PURITY 

A12 

17.2* 

7075- T6 

OLR  1.8 

18.1* 

OCR  4000 

i 

19.4* 

A13 

11.2 

7075-T7 

CA 

HI  Cu 

A14 

13.7 

7050-T7 

CA 

A16 

9.5 

7075-T7 

CA 

A17 

3.7 

7075-T6 

CA 

A20 

9.5 

7075-T7 

OLR  1.4  OCR  4000 

A21 

3.8 

7075  T6 

OLR  1.4  OCR  4000 

A26 

8.6 

7050-T7 

OLR  1.8  OCR  4000 

A27 

10.2 

7050-T7 

OLR  1.8  OCR  4000 

A30 

17.7* 

7075-T6 

OLR  1.8  OCR  4000 

A31 

0.91 

7010-T7 

OLR  1.8  OCR  8000 

HI  PURITY 

A3  2 

3.4 

7475-T7 

OLR  1.8  OCR  8000 

A33 

1.6 

7475-T7 

OLR  1.8  OCR  8000 

HI  CU 

A34 

2.75 

7010-T7 

OLR  1.8  OCR  8000 

A35 

3.20 

7050-T7 

OLR  1.8  OCR  8000 

LO  PURITY 

A36 

4.0 

7075-T7 

OLR  1.8  OCR  8000 

HI  Cu 

A37 

2.55 

7050-T7 

OLR  1.8  OCR  8000 

A39 

5.4 

7075-T7 

OLR  1.8  OCR  8000 

A40 

4.3 

7075-T6 

OLR  1.8  OCR  8000 

A41 

3.1 

7075-T6 

OLR  1.8  OCR  8000 

A42 

3.1 

7475-T7 

OLR  1.8  OCR  8000 

A43 

3.2 

7475-T7 

OLR  1.8  OCR  8000 

A44 

19.2 

7050-T7 

LHPL 

A45 

15.0 

7050-T6 

LHPL 

A46 

18.2 

7075-T7 

LHPL 

A47 

11.6 

7075-T6 

LHPL 

i 

LOAD  TYPES; 

CA  Constant  Amplitude 

OCR  Occurrence  ratio  of  1  overload  per  4000  or  8000  cycles. 

OLR  Overload  ratio  1.4,  IS 

LHPL  Lo-Hi-Lo  Programme  Loading. 

*  At  a  crack  length  2.23  cm. 


CASE: 


Lives 


TABLE  II 

EXPERIMENTAL  DESIGNS  FOR  ANALYSIS 
Rearrangement  of  data  from  Table  I 

ONE  WAY  CLASSIFICATIONS 


A7 (a)  A8 


A26 

A31 

A35 

8.6 

0.91 

3.20 

10.2 

2.75 

2.55 

Logarithmic  standard  deviation  s  *  0.1449  (6  d.f. ) 


CASE: 

Al2 (a) 

A32 

Lives 

17.2 

3.4 

cm 

18.1 

1.6 

19.4 

3.1 

17.7  * 

3.2 

s  -  0.1097  (6  d.f.) 


(a)  First  entry  in  Table  I 
*  Measured  at  a  crack  length  2.28  cm. 


TWO  WAY  CLASSIFICATIONS  (WITHOUT  INTERACTION) 


Lo-Hi-Lo  Programme  Load 


Alloy 


7050 

7075 


|  Temper 

Load  Programme 

CA 

1.4/4000 

1 

.8/8000  ] 

T6 

T7 

3.7 

9.5,  11.2 

3.8 

9.5 

a 

.3,  3. 

.0,  5. 

1 

S  =■  0.124)  (5  d.f.) 

A17,  A21,  A40,  A41,  Al6,  A  13,  A20,  A36, 
A39  by  rows. 


Temper 

T7 

T6 

19.2 

18.2 

15.0 

11.6 

s  -  0.0442  (1  d.f.) 
A44  to  A47  by  rows 


POOLING 


Degrees  of  freedom  *6  +  6  +  5  +  l«l8 

Pooled  sums  of  squared  errors  -  0.1259  6789  +  0.0721  8733  4-  0.0769  8458 

+  0,0019  5379 


Variance  estimate  ■  0.2770  9359/16  -  0.01C~  9409 
a  -  0.1240  7291  (18  d.f.) 


APPENDIX 


ANALYSIS  OF  NON-ORTHOGONAL  DATA 


■ 

!  Bi 

B2 

B3 

T6 

■ 

3.8 

BB 

T7 

BB 

9.5 

BO 

Bi  -  Constant  amplitude 
B2  -  OLR  1.4 
B3  -  OLR  1.8 


The  table  above  reproduces  the  data  in  question  from  Table  2. 

As  shown,  it  may  be  explained  by  the  model, 

■  ti  ■  bi  *  *uk 

where  the  y's  are  logarithmic  data  and  the  usual  overall  mean  has  been 
subsumed  into  the  other  parameters.  Nevertheless  there  are  still  six 
parameters  for  a  model  of  rank  five  so  that  the  normal  equations  are 
singular.  Since  the  lack  of  orthogonality  precludes  the  usual  analysis 
we  shall  solve  the  normal  equations  directly  with  allowance  for  the 
singularity.  This  amounts  to  using  the  pseudo- inverse  of  the 
coefficient  matrix. 

Each  of  the  entries  in  the  design  can  be  explained  by  one  of 
the  equations 


'■  i 

*  ? 

i 

f 

S' 


i 

* 

f 


T6  T7  B1  B2  B3 

1.1  . 

1  1 

1  .1 

1  .1 

.  1  1 

1  1 

1  .  1 

1 

.  1 


— i 

Te 

■  log 

~  3.7 

t7 

3.8 

Bl 

4.3 

B2 

3.1 

B3 

9.5 

11.2 

9.5 

1 

1 


9.5 

4.0 

5.4 


(b) 


APPENDIX  (COWTD  . ) 


or 


X6  *  y 


which  arc  simply  the  model  applied  to  each  of  the  entries  in  the  data* 


Premultiplication  by  provides  the  normal  equations 


— 

4.112 

T6 

SC 

2.2728 

.5212 

t7 

4.3391 

12  3.. 

Bl 

2.5951 

11.2. 

b2 

1.5575 

2  2  .  .  4 

b3 

2.4593 

.  -  — 

—  ~ 

L. 

Sum  check 

10.2728 

14.3391 

8.5951 

5.5575 

10.4593 


Next,  the  triangular  factors  provided  by  Crout*s  method  and  the 
corresponding  right  hand  side  combine  in  the  auxiliary  matrix 


4.0 

- 

.25 

.25 

.5 

0.5682  0387 

2.5682  0387 

. 

5.0 

.4 

.2 

.4 

0.8678  2380 

2.3678  2380 

1.0 

2.0 

1.95 

-.333* 

-  .  666 

0.1493  8045 

.1493  8045 

1.0 

1.0 

-.65 

1.333 

-1.0 

0.1639  3262 

.1639  3262 

2.0 

2.0 

-1.3 

-1.333 

0 

0/0  use  zero 

-2.9/0  X  KT1 

in  which  th^  pivot  £  ultimately  reflects  the  linear  deoendence  of  the 
model  through  the  fact  that  the  B3  column  of  X  depends  on  the  first 
four  columns. 


In  the  back  solution  the  solution  for  B3  is  arbitrary  and  for 
convenience  here  we  use  B3  *  0.  The  general  solution  may  be  found  by 
adding  another  back  solution  column  starting  from  b3  «  1.0  and  combining 
the  two  with  B3  as  an  arbitrary  coefficient  of  the  second  set  of 
solutions . 


As  in  fully  ranked  models  the  general  linear  hypothesis 
indicates  that  the  regression  sum  of  scruares  *  (X^y) ,  the  sum 


of  fitted  parameters  multiplied  by  their  respective  right  members 
in  the  normal  equations.  Here  this  is  more  conveniently  obtained  from 


APPENDIX  (CONTD.) 


the  pivots  and  right  member  of  the  auxiliary  equation  in  the  typical 
form 


RgS  =  4^.  5682)  2  +  5^.  8678) 2  +  ...0(0) 2 
*  5.1363  5845 


on  four  degrees  of  freedom,  the  rank  of  the  model.  This  of  course 
includes  the  effect  of  the  overall  mean 


C  = 


(lyijk)  /Total  number 


4.8575  196°  (1  d.f.) 


and  the  difference  between  these  numbers  are  marginal  T^  or  Bj  effects 
of  amount 


*  5.14  '  4.86  =  0.2788  3876;  (3  d.f.) 

further  reduction  needs  other  techniques^.  One  of  these  take  advantage 
of  the  elimination  inherent  in  the  auxiliary  matrix. 

From  the  last  three  {i.e.  2)  rows 

Block  SS  »  1. 95(  .1493)  2  +  1.  33*( .  1639)  2 
=  0.0793  4516?  (2  d.f.) 

The  first  two  rows  similarly  provide  and  overall  mean  scuares 
Treatment  SS  +  C  ■  5.0570  1324;  (1  d.f.) 

Consistency  with, 

^  *  Block  SS  +  Treatment  SS 


checks  the  arithmetic  of  this  stage. 


For  the  error, 


E  * 

SS 


5.2133  4303  -  4.8575  1969 
0.0769  8453?  (5  d.f.) 


corresponding  to  the  variance  estimate  0.0153  9692  or 
s  *  0.1240  8*35 


on  five  degrees  of  freedom. 


(d) 

APPENDIX  (COMTD.) 


It  is  now  convenient  to  display  these  results  in  the  ANOVA. 


Source 

Degrees  of 

Sum  of 

Mean 

F 

freedom 

squares 

square 

Temper 

1 

0.1994 

9355 

Same 

12.96 

(2.5%) 

Blocks ,  Bi 

2 

0.0793 

4516 

2.58 

(20%) 

®TB 

3 

0.2738 

3876 

Mean  C 

1 

4.8575 

1969 

Error 

5 

0.0769 

3458 

0.0153  9692 

TOTAL 

9 

5.2133 

4303 
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